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I. INTRODUCTION 

In a recent work [l[ an alternative gravity theory with 
nonminimal coupling has been proposed. The theory ex- 
hibits - in contrast to Einstein's theory of gravitation - a 
direct coupling between the matter fields and the curva- 
ture scalar on the Lagrangian level. In contrast to other 
theories with extended dynamics in the gravitational sec- 
tor, such a nonminimal coupling prescription leads to a 
modification of the equations of motion. The latter were 
analyzed for specific choices of the matter Lagrangian, 
e.g., a perfect fluid 0- 

Here we present a systematic derivation of the equa- 
tions of motion for arbitrarily structured test bodies. The 
method employed by us is not tied to a specific choice of 
the matter Lagrangian and therefore generalizes previ- 
ous findings. The propagation equations for pole-dipole 
particles are worked out with the help of a multipole 
method. In particular it allows us to provide the form 
of the extra force terms - "extra" in comparison to the 
case with minimal coupling - entering the equations of 
motion. Therefore, our findings should be taken into ac- 
count in the context of the recent controversy, see 0, Q , 
regarding the appearance of such an extra term. Fur- 
thermore, our results allow for direct a comparison with 
the motion of test bodies in other gravity theories, in 
particular the coupling between material and geometri- 
cal quantities becomes evident. 



II. THE MODEL UNDER CONSIDERATION 

In [l[ an extended version of a so-called f{R) gravity 
theory was considered - for earlier works on this subject 
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see also [j, l5| , as well as [Q| for an extension which covers 
bodies with spin/polarization. Gravity theories in which 
the usual Einstein-Hilbert Lagrangian is replaced by an 
arbitrary function of the curvature scalar have attracted 
a lot of attention during the last few years see, e.g., the 
reviews 0, [1| and references therein. The /(i?)-scenario 
is generalized even further in [l| by the introduction of 
a nonminimal coupling term on the Lagrangian level. In 
particular the following Lagrangian was put forward 

itot = ^/l (i?) + [1 + A/2 (i?)] L^af (1) 

Here /i and /2 are arbitrary functions of the curvature 
scalar R and Lmat is the matter Lagrangian. The non- 
minimal coupling of matter and gravity is controlled by 
the constant A. The general field equations - in terms 
of the functions /i and /2 and their derivatives - are 
given in [1,] ; their explicit form is irrelevant for the subse- 
quent analysis though. Theories of the above mentioned 
kind, in particular with a nonminimal coupling between 
the curvature scalar and a scalar field, have also been 
considered before in a cosmological context see, e.g., the 
review 

In contrast to standard f{R) gravity theories, the last 
term in ([1]) leads to a modification of the equations of 
motion. As was already shown in (5) of [1] the usual 
conservation law - as, for example, found in general rel- 
ativity - is replaced by 

t + A/2 

Here F2 (i?) := (i/2 (R) /dR denotes a shortcut for 
derivatives of the unspecified function /2 (i?) of the 
curvature scalar and the energy-momentum tensor of 
matter is defined in a standard way by y^—gTij :— 

Of course from ^ it becomes immediately apparent, 
that the equations of motion of the theory under con- 
sideration differ from the ones of general relativity if the 
right-hand side (rhs) of ^ is nonvanishing. In \V\ it is 
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suggested that a deviation of this kind due to the non- 
minimal couphng term may play a role for the observed 
flatness of rotation curves IHaIiII or in the context of the 
so-called Pioneer anomaly [12]. 

In their study of the motion of test bodies, the au- 
thors of [H made an explicit assumption for the energy- 
momentum tensor of matter entering ([2|) , which they as- 
sumed to be of perfect fluid form. Also in [isj a similar 
analysis is carried out which is based on a specific choice 
for the energy-momentum tensor of the system. 

Here we extend the previous analysis by deriving the 
general form of the equations of motion for test bodies. 
Our analysis, which relies on a well established multi- 
pole approximation technique, is independent of a spe- 
cific choice of energy-momentum tensor. 

Without going into historical detail we only mention 
that in the context of general relativity the method 
- and variations of it - was utilized in the works 

[3, Hi, M, [13, M, M, H ISil, [13, M, H M, 111], it 

has also been successfully applied in alternative gr avity 
theories, see [13, HE [23 and more recently in [sO, [3l| . 
Note that reference [30| also contains a short timeline 
of works. The model under consideration in the present 
work does not belong to the very general class of grav- 
itational models analyzed in (30j due to its nonminimal 
coupling prescription. 



Here Sx"^ := a;" — and Y{t) parametrizes the world 
line of the body. Note that here, and in the following, 
only the last two indices of the integrated quantities are 
spacetime indices, and the bi ■ ■ ■ bn label the multipole 
order of a current. The integrals in ([5]) are taken over a 
3-dimensional slice at a time t, over the world tube 
of the test body. We use the condensed notation 



/ fix)d'x. 



(6) 



With the definitions in ([5]) the integrated version of the 
conservation law ([3]) takes the following general form: 



dt 



/3=1 



here an inverted circumflex, e.g. "6^" , indicates the omis- 
sion of an index from a list and v"" := dV^/dt 



IV. PROPAGATION EQUATIONS FOR 
POLE-DIPOLE TEST BODIES 



III. INTEGRATED CONSERVATION LAW 

The multipole scheme employed by us relies on the 
integration of the conservation law over the world tube 
of the test body, see also HO] for details. To begin with, 
we rewrite ([2]) as follows: 



In this section we work out the equations of motion 
for pole-dipole test bodies. For such bodies only the mo- 
ments T*"' , t''*"' , S*"* , and s'^*"' are nonvanishing. With the 
expansion of geometrical quantities around the worldline 
Y{t) of the test particle into a power series in Sx"' 



ij -'-'mat 



(3) 



Here we introduced a scalar function A (i?) :— 
log [1 -|- A/2 (R)]- In the following we are going to denote 
derivatives of this function simply by Ai := diA,Aij := 
dfjA, etc. Raising the indices and rewriting the covariant 
derivative in ([3]) yields 



(4) 



In the last equation we introduced the quantity S*-' :— 
^'■'imat as a shortcut. Densities of different quantities 
are denoted by a tilde " ". We define the integrated 
multipole moments of the matter quantities S'^ and T'^ 
as follows: 



T 



—bi---b„ij 



(5) 



R\ 



R\y+Sx'' Ra 

1 

+ 2 



Sx^'Sx" R^ab\Y 



= r 



«j \y 



+ Sx" r 



ij ,a|y 



(8) 



the integrated conservation law ([7]) yields the following 
set of propagation equations: 



d_—m 
dt 



d_—aiO 

dt 



^ —cib 



—cd 



i cd .h-L 



= 1 ~ V 1 -icrfJ 



—aib TT^aib\ , 

T ]Ab, 



= 



(9) 

(10) 
(11) 



Here we suppressed the dependencies on the points at 
which the quantities are evaluated. 
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V. PROPAGATION EQUATIONS REWRITTEN 

To allow for a better comparison with the result in [13] , 
we bring - (fTT|) into a form which closely resembles 
the form of the covariant equations (5.3) and (5.7) in 

We start with the following redefinitions of the inte- 
grated moments: 

Ar" u'^T', M-"' := -u^'T"'. (12) 

Here we introduced m° :— dV^/ds for the velocity and 
the parameter s denotes proper time. Furthermore one 
should note that M°°'' = due to our choice of the inte- 
gration domain over hypersurfaces with t = const. Anal- 
ogously to [17i] we introduce the spin as follows: 

5"'':=r'°-T'"°. (13) 
From these definitions we can immediately infer that 

With the propagation equation (jlip becomes 

u° (iVP" + AP"^) = u'AP''^ + u^M"'^. (15) 

Cyclic permutation of the indices in equation (|15p and 
subtraction of the second from the sum of the first and 
the third of the permutations yields 



2M''^J ^ - (m^S-" + u^S'^) + -Q {S^^u" + S^^u^) . (16) 
With the definition 



iV-" := (17) 



the second propagation equation (jlOp becomes 



M"'^ d f AP' 



vP ds \ 
_ (TV^'^f' + M'''^^) Ab. (18) 

If we introduce the following "generalized momentum" 

J_ (^jao ^ r,/w'*5=°) , (19) 
equation (fT^ can be recast into 

_ (iV^'^b + M"'') Af,. (20) 
Here we used the definition of the covariant derivative 

n qab 1 cab 

■=^ + Tcd'^S^V + TjS'^^u^. (21) 
Us ds 

Taking the skew symmetric part of (I20|) yields 

1 DS" 



2 Ds 



0, 



(22) 



and the derivative with the hat is defined as follows 
DS'" DS"" 



Ds 



Ds 



+ ( TV'*'"' - TV*"" + M° 



NP' 



A. (23) 



If we contract ((22t with Ui and make the same choice for 
the mass as in [17], namely 



m := fj, Ui, 



we obtain 



mu + Ui- 



DS" 



Ds 



(24) 



(25) 



Substituting this equation back into yields a very 
compact version of the second propagation equation, 
namely 



DS'"" , DS""" „ DS"" „ 
-f u Uc—p: u u^— — = 0, 



Lis 



Ds 



Ds 



which should be compared to (5.3) in [I] 
Finally, with the definition 



N"'' := u°E 
the first propagation equation ^ becomes 



(26) 



(27) 



M 



aO 



ds \ 



a n ,fbcd 
i cd .b^d 



cab 



AP"''] A 



be- 



(28) 



With the help of ([^0]) and and by using the deriva- 
tives as defined in (pij) and (^5]) we can bring (^5)1 into 
its final form 



D 
ITs 



DS"" 

nu + Ur 

Ds 

( Afcdfc _|_ ]\.[cdb\ 



+ ^S^'^u'^Rbcd" = {N"'' - NP^) Ai 



{S^Akc + Fe/Afc) . (29) 
This result should be compared to equation (5.7) in pTzj . 

VI. PHYSICAL CONSEQUENCES 

Our main result is embodied in the propagation equa- 
tions as given in (pS)) and ([211). These equations should 
be compared to the well-known ones for pole-dipole test 
par ticles in general relativity as given in (5.3) and (5.7) of 
TT] , nowadays usually called the Mathisson-Papapetrou- 
Dixon (MFD) equations. 

Our result in the case of the nonstandard gravity the- 
ory clearly shows the additional terms which arise due 
to the nonminimal coupling. The structure of the gen- 
eralized propagation equation for the spin, i.e. equation 
(|26)) . is very similar to the classic result in (5.3) of p7j . 
In (pS)) the new derivative from (P^)) takes into account 
the extra terms which arise due to the nonminimal cou- 
pling prescription on the Lagrangian level. The second 



4 



of the MPD equations can be easily recovered from (l26l) 
by replacing the "hatted" derivative by the standard one 
given in ((2T|) . Also the first generalized propagation (|29|) 
allows for a quick recovery of the first MPD equation for 
the momentum. As soon as one switches off the nonmin- 
imal coupling, the rhs of ((29|) vanishes, and the hatted 
derivative on the left-hand side (Ihs) is replaced by the 
standard derivative. 

It is very interesting to note that, even when we confine 
ourselves to nonspinning particles - i.e. test particles for 
which M°-'"^ vanishes ~ the rhs of is nontrivial and 
the generalized momentum /i*, which enters on the Ihs of 
([29)) . is given by 

= mu' + {N^''' - N'^^) UcAb, (30) 

yielding fj} ^ mu^ . 

Furthermore, if we do not allow for any kind of dipole 
contribution, neither via TV''^ nor M"'"^ - i.e. we con- 
sider only pole particles, the first propagation equation 
becomes 

^ {mu') - (7V'^ - AP'') At. (31) 
Us ' ^ ' 

Hence, we still have a nonstandard contribution on the 
rhs of the first equation of motion, with a direct coupling 
between the monopole moments of the matter currents 
and the background geometry as described by the deriva- 
tive of the function A(TC). 

In other words, even single-pole test particles do not 
move along geodesies in the theory under consideration. 
Of course this is in contrast to the standard result in the 
theory of general relativity. 

One should note that in the case of single-pole particles 
- and by using and (pS)) - equation (pij) can be 
brought into the following form: 

{tuu^ (1 -f A/2)] = TV^" (1 + A/2)^, . (32) 

From this equation one can immediately read off the 
additional contribution due to the nonminimal coupling 
terms as embodied by f2{R)- On the Ihs of the single- 
pole equation of motion, these terms induce an "effec- 
tive mass" . On the rhs of p2|) , the nonminimal coupling 
procedure ~ remembering TV"'' = g'^''N - leads to some 
"effective pressure" term. 

Equation p2p provides an interestin g in terpretation of 
the aforementioned Pioneer anomaly [l2] as a result of 
the nonminimal coupling of matter and gravity. From 
([32|) it is obvious that no extra forces act on a test body 
for gravitational field configurations with vanishing cur- 
vature scalar R = 0. 

However, this picture changes as soon as one considers 
spacetimes with nonvanishing curvature scalar, for ex- 
ample as encountered in cosmological solutions. Albeit 
not an exact solution in the context of the theory un- 
der consideration, the well-known Friedmann-Lemaitre- 
Robertson- Walker (FLRW) metric may serve as an il- 
lustrating example. For a linear choice of the function 



/2 - or, equivalently, for the first term in a series ex- 
pansion of /2 - the extra force entering ([5^ would be 
directly proportional to the curvature scalar. Hence, in a 
FLRW background one would obtain an additional con- 
tribution due to the nonminimal coupling prescription 
which is proportional to 

f2 ^ Rflrw = 6 (^H^ {I -q) + ^y (33) 

Here we made use of the standard definition of the 
Hubble rate H := a/ a and the deceleration parameter 
q := —aa/a?. The scale factor is a function of time only, 
i.e. a — a{t), and the constant k determines the spatial 
curvature in the FLRW metric. 

Hence, if we consider the motion of a test body, e.g. 
of the Pioneer spacecraft, taking into account the cos- 
mological expansion of the spacetime background, then 
one would obtain a correction to its acceleration - com- 
pared to geodesic motion of general relativity - which 
is proportional to the Hubble rate H and its derivatives 
(which can probably be neglected). This is an interest- 
ing result which qualitatively agrees with the measured 
anomalous acceleration of the Pioneer spacecraft which 
is approximately cHq ^ 7 x 10"-'^'' ms~^. A quantitative 
analysis could further provide estimates for the nonmini- 
mal coupling constant A and for the form of the function 

VII. CONCLUSION 

We worked out the explicit form of the equations of 
motion for a gravitational theory ([T]) with nonminimal 
coupling. Our results extend previous works on the equa- 
tions of motion of such theories and are independent of 
the specific form of the energy-momentum tensor. 

The coupling between geometric and matter quantities 
becomes apparent in our framework and should be taken 
into account in the systematic testing of the theory. Fur- 
thermore, our analysis confirms ~ in a very general way 
- the nongeodesic motion of single-pole test bodies. 

The analysis in this work also applies to several other 
models with nonminimal coupling. A direct compari- 
son can be made by a simple remapping of the function 
A{R) in the present analysis. Our results also allow for a 
straightforward comparison with the equations of motion 
in other alternative gravity theories [23, [H, |2^ [s^, HH , 
which do not belong to the class of f{R) gravity theories. 

Furthermore, we have shown that within a nonminimal 
coupling scheme one can expect an additional accelera- 
tion of bodies due to the global influence of an expanding 
universe. This is particularly interesting in the context 
of the observed anomalous acceleration of the Pioneer 
spacecraft, since our results indicate that the "cosmolog- 
ical order" of this effect is more than a mere coincidence, 
but could be ascribed to the nonminimal coupling. 

It would be interesting to carry out a systematic study 
of the motion of test bodies in specific background space- 
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times. Such an analysis should lead to very tight con- 
straints of the free parameters of the nonminimal cou- 
pling model. 
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